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Astrophysical systems known as inspiralling compact binaries are among the most interesting sources to hunt for 
gravitational radiation in the future network of laser-interferometric detectors, composed of the large-scale interfer- 
" ometers VIRGO and LIGO, and the medium-scale ones GEO and TAMA (see the books @-| for reviews, and the 
contribution of B Schutz in this volume). These systems are composed of two compact objects, i.e. gravitationally- 
condensed neutron stars or black holes, whose orbit follows an inward spiral, with decreasing orbital radius r and 
increasing orbital frequency uj. The inspiral is driven by the loss of energy associated with the gravitational- wave 
emission. Because the dynamics of a binary is essentially aspherical, inspiralling compact binaries are strong emitters 
of gravitational radiation. Tidal interactions between the compact objects are expected to play a little role during 
most of the inspiral phase; the mass transfer (in the case of neutron stars) does not occur until very late, near the 
final coalescence. Inspiralling compact binaries are very clean systems, essentially dominated by gravitational forces. 
Therefore, the relevant model for describing the inspiral phase consists of two point-masses moving under their mu- 
tual gravitational attraction. As a simplification for the theoretical analysis, the orbit of inspiralling binaries can be 
t-H i considered to be circular, apart from the gradual inspiral, with a good approximation. At some point in the evolution, 
there will be a transition from the adiabatic inspiral to the plunge of the two objects followed by the collision and final 
merger. Evidently the model of point-masses breaks down at this point, and is to be replaced by a fully relativistic 
numerical computation of the plunge and merger (see the contribution of E Seidel in this volume) . 

Currently the theoretical prediction from general relativity for the gravitational waves emitted during the inspiral 
phase is determined using the post-Newtonian approximation (see [|]|5| for reviews). This is possible because the 
t-H , dynamics of inspiralling compact binaries, though very relativistic, is not fully relativistic: the orbital velocity v is 
always less than one third of c (say). However, because 1/3 is far from negligible as compared to 1, the gravitational- 
j radiation waveform should be predicted up to a high post-Newtonian order. In particular, the radiation reaction 
Q-i onto the orbit, which triggers the inspiral, is to be determined with the maximal precision, corresponding to at least 
the second and maybe the third post-Newtonian (3PN, or 1/c 6 ) order Notice that the zeroth order in this 

(^j), post-Newtonian counting corresponds to the dominant radiation reaction force (already of the order of 2.5PN relative 
to the Newtonian force), which is due to the change in the quadrupole moment of the source. Actually, the method 
.^h ! is not to compute directly the radiation reaction force but to determine the inspiral rate from the energy balance 
equation relating the mechanical loss of energy in the binary's centre of mass to the total emitted flux at infinity. 

The implemented strategy is to develop a formalism for the emission and propagation of gravitational waves from 
a general isolated system, and only then, once some general formulae valid to some prescribed post-Newtonian order 
are in our hands, to apply the formalism to compact binaries. Hence, we consider in this paper a particular formalism 
applicable to a general description of matter, under the tenet of validity of the post-Newtonian expansion, namely 
that the matter should be slowly moving, weakly stressed and self-gravitating. Within this formalism we compute 
the retarded far field of the source by means of a formal post-Minkowskian expansion, valid in the exterior of the 
source, and parametrized by some appropriately defined multipole moments describing the source. From the post- 
Minkowskian expansion we obtain a relation (correct up to the prescribed post-Newtonian order) between the radiative 
multipole moments parametrizing the metric field at infinity, and the source multipole moments. On the other hand, 
the source multipole moments are obtained as some specific integrals extending over the distribution of matter fields 
in the source and the contribution of the gravitational field itself. The source moments are computed separately up to 
the same post- Newtonian order. The latter formalism has been developed by Blanchet, Damour and Iyer JH-H. More 



recently, a different formalism has been proposed and implemented by Will and Wiseman |L5j (see also |16 171]). The 



two formalisms are equivalent at the most general level, but the details of the computations are quite far apart. In the 
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second stage, one applies the formalism to a system of point-particles (modelling compact objects) by substituting 
for the matter stress-energy tensor that expression, involving delta-functions, which is appropriate for point-particles. 
This entails some divergencies due to the infinite self-field of point-particles. Our present method is to cure them 
systematically by means of a variant of the Hadamard regularization (based on the concept of 'partie finie') Jis|jl9t| . 

In this paper, we first analyse the binary inspiral gravitational waveform at the simplest Newtonian approximation. 
Notably, we spend some time describing the relative orientation of the binary with respect to the detector. Then we 
compute, still at the 'Newtonian order' (corresponding, in fact, to the quadrupole approximation), the evolution in the 
course of time of the orbital phase of the binary, which is a crucial quantity to predict. Next, we review the main steps 
of our general wave-generation formalism, with emphasis on the definition of the various types of multipole moments 
which are involved. At last, we present the result for the binary inspiral waveform whose current post-Newtonian 
precision is 2PN in the wave amplitude and 2.5PN in the orbital phase (that is 1/c 5 beyond the quadrupole radiation 
reaction). However, since our ultimate aim is to construct accurate templates to be used in the data analysis of 
detectors, it is appropriate to warm up with a short review of the optimal filtering technique which will be used for 
hunting the inspiral binary waveform (see pcfl for an extended review). 



II. SUMMARY OF OPTIMAL SIGNAL FILTERING 



Let o(t) be the raw output of the detector, which is made of the superposition of the useful gravitational-wave 
signal hit) and of noise n(t): 

o(t) = h(t) + n(t). (2.1) 
The noise is assumed to be a stationary Gaussian random variable, with zero expectation value, 

n(t) = 0, (2.2) 
and with (supposedly known) frequency-dependent power spectral density S n (u>) satisfying 



n(u)n*(u f ) = 2it6(uj - u')S n (u>), (2.3) 



where h(u>) is the Fourier transform of n(t). In ( |2.2| ) and Q2.31) , we denote by an upper bar the average over many 
realizations of noise in a large ensemble of detectors. From (|2.3j), we have S n (uj) — S*(tu) = S„(—lo) > 0. 

Looking for the signal h(t) in the output of the detector o(t), the experimenters construct the correlation c[t) 
between o(t) and a filter q(t), i.e. 

/+oo 
dt'o(t')q(t + t'), (2.4) 
-oo 

and divide c(t) by the square root of its variance (or correlation noise). Thus, the experimenters consider the ratio 

c(t) f + °° ^5(uj)q*(Lu)e iujt 
°wK*) = — — — 2 : z = ~ n72> ( 2 - 5 ) 



(c 2 (*)-c(*))V 2 (f+~^S„ Hk >)| 



where o(u>) and q(u>) are the Fourier transforms of o(t) and q(t). The expectation value (or ensemble average) of this 
ratio defines the filtered signal-to-noise ratio 

( + °° p-h{w)q*(u)e lult 

Plq](t) - <r[q](t) = , °° 27r /2 - (2.6) 



The optimal filter (or Wiener filter) which maximizes the signal-to- noise (2.6) at a particular instant t = (say), is 
given by the matched filtering theorem as 

(2-7) 

b n {u>) 
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where 7 is an arbitrary real constant. The optimal filter (2.7) is matched on the expected signal h(uj) itself, and 
weighted by the i nver se of the power spectral density of the noise. The maximum signal to noise, corresponding to 
the optimal filter (2.7), is given by 



P = 



dcu \h(tu)\' 
2tt S n (u;) 



2\ 1/2 



(h,h) 



1/2 



(2.8) 



This is the best achievable signal-to-noise ratio with a linear filter. In (2.S), we have used, for any two real functions 
f(t) and g(t), the notation 



(f,g) 



+00 



da; f(uj)g*{uj) 
2-7T S n {uo) 



(2.9) 



for an inner scalar product satisfying (/,<?) = (/,<?)* = (<?,/}• 

In practice, the signal h(t) or h{uj) is of known form (given, for instance, by (3.9)-( j3T4| ) later) but depends on an 
unknown set of parameters which describe the sour ce o f radiation, and are to be measured. The experimenters must 
therefore use a whole family of filters analogous to (2.7) but in which the signal is parametrized by a whole family of 
'test' parameters which are a priori different from the actual source parameters. Thus, one will have to d efine and 
use a lattice of filters in the parameter space. The set of parameters maximizing the signal to noise ( |2.6| ) is equal, 
by the matched filte ring theorem, to the set of source p aram eters. However, in a single detector, the experimenters 
maximize the ratio (2.5) rather than the signal to noise (2.6), and therefore make errors on the determination of the 
parameters, depending on a particular realization of noise in the detector. If the signal-to-noise ratio is high enough, 
the measured values of the parameters are Gaussian distributed around the source parameters, with variances and 
correlation coefficients given by the covariance matrix, the computation of which we now recall. Since the optimal filter 
(2.7) is defined up to an arbitrary multiplicative constant, it is convenient to treat separately a constant amplitude 
parameter in front of the signal (involving, in general, the distance of the source). We shall thus write the signal in 
the form 



h(cj;A, X a ) = A k(uj;X a ), 



(2.10) 



where A denotes some amplitude parameter. The function k depends only on the other parameters, collectively 
denoted by A a where the label a ranges on the values 1,...,N. The family of matched filters (or 'templates') we 
consider is defined by 



q{uJ\tK) = 7 



, fc(u; t A a ) 
S n (u}) 



(2.11) 



wher e t A a is a set of test parameters, assumed to be al l in dependent, and 7' is arbitrary. By substituting (2.11) into 
(2.5) and choosing t = 0, we get, with the notation of (|2.9|), 



(o,k( t X)) 



(fc(*A),fc( t A))i/2 



(2.12) 



(Note that a is in fact a function of both the parameters A a and tA a .) Now the experimenters choose as their 
best estimate of the source par ameters A a the measured parameters m X a which among all the test parameters tX a 
(independently) maximize (ET2j), i.e. which satisfy 



da 

d t X a 



( m X) =0, a = 1, . 



,N. 



(2.13) 



Assuming that the signal to noise is high enough, we can work out (2.13) up to the first order in the difference between 
the actual source parameters and the measured ones, 



As a result, we obtain 



5X a = X a 



dh 



5X a = C ab { -(n, tpt) + 



dX (h, h) ' dXb' 



(2.14) 



(2.15) 
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where a summation is understood on the dummy label b, and where the matrix C a b (with a, b 
inverse of the Fisher information matrix 



1, . . . , N) is the 



dh dh 
dX a ' dX b 



1 



(h, h) \ ' d\ a 



(2.16) 



(we have C a bT>b c = S ac ). On the right-hand sides of ( 2.1 5| ) and flS.lq 1 ), the signal is equally ( with t his approximation) 
parametrized by the measured or actual parameters. Since the noise is Gaussian, so are, by ( 2.1 5| ) , the variables <5A a 
(indeed, 6X a result from a linear operation on the noise variable). The expectation value and quadratic moments of 
the distribution of these variables are rea dily o btai ned from the facts that (n, f)—0 and (n, f) (n, g) = (f, g) for any 
deterministic functions / and g (see (2.2) and ( |2.3|) ). We then obtain 



SXa 

5X a SXb 



0. 

Cab- 



(2.17) 



Thus, the matrix C a b (the inverse of ( p. 16 )) is the matrix of variances and correlation coefficients, or covariance matrix, 
of the variables <5A a . The probability distribution of SX a reads as 



P(SX a ) = 



v /(27r) JV + 1 detC 



CXp { -\v ab 5X a 5Xi 



(2.18) 



where detC is the determinant of C a b- A similar analysis can be done for the measurement of the amplitude parameter 
A of the signal. 



III. NEWTONIAN BINARY POLARIZATION WAVEFORMS 



The source of gravitational waves is a binary system made of two point-masses moving on a circular orbit. We 
assume that the masses do not possess any intrinsic spins, so that the motion of the binary takes place in a plane. 
To simplify the presentation we suppose that the centre of mass of the binary is at rest with respect to the detector. 
The detector is a large-scale laser- interferometric detector like VIRGO or LIGO, with two perpendicular arms (with 
length 3 km in the case of VIRGO). The two laser beams inside the arms are separated by the beam-splitter which 
defines the central point of the interferometer. We introduce an orthonormal right-handed triad (X, Y, Z) linked with 
the detector, with X and Y pointing along the two arms of the interferometer, and Z pointing toward the zenithal 
direction. We denote by n the direction of the detector as seen from the source, that is, — n is defined as the unit 
vector pointing from the centre of the interferometer to the binary's centre of mass. We introduce some spherical 
angles a and f3 such that 

— n = X sin a cos /3 + Y sin a sin (3 + Z cos a. (3-1) 

Thus, the plane (3 = constant defines the plane which is vertical, as seen from the detector, and which contains the 
source. Next, we introduce an orthonormal right-handed triad (x, y, z) which is linked to the binary's orbit, with 
x and y located in the orbital plane, and z along the normal to the orbital plane. The vector x is chosen to be 
perpendicular to n; thus, n is within the plane formed by y and z. The orientation of this triad is 'right-hand' with 
respect to the sense of motion. We denote by i the inclination angle, namely the angle between the direction of the 
source or line-of-sight n and the normal z to the orbital plane. Since z is right-handed with respect to the sense of 
motion we have < i < n. Furthermore, we define two unit vectors p and q, called the polarization vectors, in the 
plane orthogonal to n (or plane of the sky). We choose p = x and define q in such a way that the triad (n, p, q) is 
right-handed; thus 

n = y sin i + z cos i, (3-2) 
p = x, (3.3) 
q = y cos i — z sin i. (3-4) 

Notice that the direction p = x is one of the 'ascending node' N of the binary, namely the point at which the bodies 
cross the plane of the sky moving toward the detector. Thus, the polarization vectors p and q lie, respectively, along 
the major and minor axis of the projection onto the plane of the sky of the (circular) orbit, with p pointing toward N 
using the standard practice of celestial mechanics. Finally, let us denote by £ the polarization angle between p and 
the vertical plane — constant; that is, £ is the angle between the vertical and the direction of the node N. We have 
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n = —X sin a cos (3 — Y sin a sin (3 — Z cos a, (3-5) 
p = X(cos £ cos a cos (3 + sin £ sin /3) 

+Y(cos £ cos a sin (3 — sin £ cos (3) — 7i cos £ sin a, (3.6) 
q = X(— sin £ cos a cos (3 + cos£ sin/3) 

+Y(— sin £ cos a sin /3 — cos £ cos /3) + Z sin £ sin a. (3.7) 

Defining all these angles, the relative orientation of the binary with respect to the interferometric detector is entirely 



determined. Indeed using (3.4) and (3.7) one relates the triad (x,y,z) associated with the source to the triad (X,Y,Z) 
linked with the detector. 

The gravitational wave as it propagates through the detector in the wave zone of the source is described by the so- 
called transverse and traceless (TT) asymptotic waveform hj^ — (g^ — (5y) TT , where g^ denotes the spatial covariant 
metric in a coordinate system adapted to the wave zone, and Sij is the Kronecker metric. Neglecting terms dying out 
like 1/R 2 in the distance to the source, the two polarization states of the wave, customarily denoted /i+ and h x , are 
given by 

h + = \(piPj - qiqj)hj?, (3.8) 
hx = +Pjqi)hJ/ : , (3.9) 

where pi and qi are the components of the polarization vectors. The detector is directly sensitive to a linear combination 
of the polarization waveforms h + and h x given by 

h(t) =F+h+{t)+Fxhx{t), (3.10) 

where F + and F x are the so-called beam-pattern functions of the detector, which are some given functions (for a given 
type of detector) of the direction of the source a, 8 and of the polarization angle £. This h(t) is th e gravitational-wave 



signal looked for in the data analysis of section [II| and used to construct the optimal filter ( 2.10 ). In the case of the 
laser-interferometric detector we have 

F + = |(1 + cos 2 a) cos 2(3 cos 2£ + cos a sin 2(3 sin 2£, (3-H) 
Fx = -§(1 + cos 2 a) cos 2/3 sin 2£ + cos a sin 2/3 cos 2£. (3.12) 

The orbital plane and the direction of the node N are fixed so the polarization angle £ is constant (in the case 
of spinning particles, the orbital plane precesses around the direction of the total angular momentum, and angle £ 
varies). Thus, the gravitational wave h(t) depends on time only through the two polarization waveforms h + (t) and 
h x (t). In turn, these waveforms depend on time through the binary's orbital phase 4>{t) and the orbital frequency 
oj(t) = d(j)(t)/dt. The orbital phase is defined as the angle, oriented in the sense of motion, between the ascending node 
iV and the direction of one of the particles, conventionally particle 1 (thus <j) = modulo 2tt when the two particles 
lie along p, with particle 1 at the ascending node). In the absence of any radiation reaction, the orbital frequency 
would be constant, and so the phase would evolve linearly with time. Because of the radiation reaction forces, the 
actual variation of <f>(t) is nonlinear, and the orbit spirals in and shrinks to zero-size to account, via the Kepler third 
law, for the gravitational-radiation energy loss. The main problem of the construction of accurate templates for the 
detection of inspiralling compact binaries is the prediction of the time variation of the phase 4>{t). Indeed, because 
of the accumulation of cycles, most of the accessible information allowing accurate measurements of the binary's 
intrinsic parameters (such as the two masses) is contained within the phase, and rather less accurate information is 
available in the wave amplitude itself. For instance, the relative precision in the determination of the distance R to 
the source, which affects the wave amplitude, is less than for the masses, which strongly affect the phase evolution 
Hence, we can often neglect the higher-order contributions to the amplitude, which means retaining only the 
dominant harmonics in the waveform, which corresponds to a frequency at twice the orbital frequency. 

Once the functions <p(t) and u>(t) are known they must be inserted into the polarization waveforms computed by 
means of some wave-generation formalism. For instance, using the quadrupole formalism, which neglects all the 
harmonics but the dominant one, we find 

h+ = T^[-^r-) (1 + cos 2 i) cos 20, (3.13) 



'R\c 3 J 

2Gp, /Gmuj\ 
1?R I c 3 J 



2/3 



(2cosz)sin2</> (3.14) 
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where R denotes the absolute luminosity distance of the binary's centre of mass; the mass parameters are given by 

m = mi + m.2; /i = ; u — — . (3.15) 

m m 

This last parameter u, introduced for later convenience, is the ratio between the reduced mass and the total mass, 
and is such that < v < 1/4 with v — > in the test-mass limit and v = 1/4 in the case of two equal masses. 



IV. NEWTONIAN ORBITAL PHASE EVOLUTION 

Let yi(t) and ya(i) be the two trajectories of the masses mi and m.2, and y = yi — y2 be their relative position, 
and denote r = |y|. The velocities are vi(t) = dyi/di, V2(i) = dy2/di and v(i) = dy/di. The Newtonian equations 
of motion read as 

dvi _ ^ m2 y. dv 2 _ ^ mi y (4 1) 

dt r 3 ' dt r 3 

The difference between these two equations yields the relative acceleration, 

dv Gm 

We place ourselves into the Newtonian centre-of-mass frame defined by 

miyi + m 2 y 2 = 0, (4.3) 
in which frame the individual trajectories yi and y2 are related to the relative one y by 



The velocities are given similarly by 



m 2 mi 

yi = — y; y2 = y. (4.4) 

m m 



m 2 mi 

vi = v; v 2 = v. (4.5) 

m m 



In principle, the binary's phase evolution cj)(t) should be determined from a knowledge of the radiation reaction 
forces acting locally on the orbit. At the Newtonian order, this means considering the 'Newtonian' radiation reaction 
force, which is known to contribute to the total acceleration only at the 2.5PN level, i.e. 1/c 5 smaller than the 
Newtonian acceleration (where 5 = 2s + 1, with s = 2 the helicity of the graviton). A simpler computation of the 
phase is to deduce it from the energy balance equation between the loss of centre-of-mass energy and the total flux 
emitted at infinity in the form of waves. In the case of circular orbits one needs only to find the decrease of the 
orbital separation r and for that purpose the balance of energy is sufficient. Relying on an energy balance equation is 
the method we follow for computing the phase of inspiralling binaries in higher post-Newtonian approximations (see 



section |V|). Thus, we write 

dE 



dt ( 46 ) 
where E is the centre-of-mass energy, given at the Newtonian order by 

£=-^P, (4.7) 

and where C denotes the total energy flux (or gravitational 'luminosity'), deduced to the Newtonian order from the 
quadrupole formula of Einstein: 

^d^d 3 ^ 

5C 5 d^ d^ ■ 1 8j 

The quadrupole moment is merely the Newtonian (trace-free) quadrupole of the source, which reads in the case of 
the point-particle binary as 



G 



(4.9) 



In the mass-centred frame (4.3) we get 



Qij = I'UIH 1 - U ij r 2 ). 



(4.10) 



The third time derivative of Qij needed in the quadrupole formula (4.8) is easily obtained. When an acceleration is 
generated we replace it by the Newtonian equation of motion (4.2). In the case of a circular orbit we get 



d 3 Qij 
dt 3 



-4— ^(yV + y j v l ) 



(this is automatically trace- free because y • v = 0). Replacing ( 4.11 ) into (O) leads to the 'Newtonian' flux 



C = 



32 G 3 m 2 [i 2 ^ 
~5 



(4.11) 



(4.12) 



A better way to express the flux is in terms of some dimensionless quantities, namely the mass ratio v given in (3.15), 
and a very convenient post-Newtonian parameter defined from the orbital frequency lo by 



/ Gmuj\ 



2/3 



(4.13) 



Notice that x is of formal order 0(1/ c ) in the post-Newtonian expansion. Thanks to the Kepler law Gm = r 3 u we 
transform (4.12) and arrive at 



32 c° 2 5 
L = ^G VX - 



In this form the only factor having a dimension is 



3.63 x 10 52 W, 



G 



(4.14) 



(4.15) 



which is the Planck unit of a power, which turns out to be independent of the Pla nck constant. (Notice that instead 
of c 5 /G the inverse ratio G/c 5 appears as a factor in the quadrupole formula (4.8).) On the other hand, we find that 
E reads simply 



-■^/ic 2 x. 



(4.16) 



Next we replace ( 4.14 ) and ( 4.16 ) into the balance equation (|4.6| ), and find in this way an ordinary differential equation 
which is easily integrated for the unknown x. We introduce for later convenience the dimensionless time variable 



5Gm 



(t c -t), 



where t c is a constant of integration. Then the solution reads 

x(t) = ir- 1/4 . 



(4.17) 



(4.18) 



It is clear that t c represents the instant of coalescence, at which (by definition) the orbital frequency diverges to 
infinity. Then a further integration yields <f> — Jwdi = — ^ J x 2 ! 3 dr, and we get the looked for result 



m = -r 

v 



5/8 



(4.19) 



where <fi c denotes the constant phase at the instant of coalescence. It is often useful to consider the number Af of 
gravitational- wave cycles which are left until the final coalescence starting from some frequency lo: 



1 

32tti/ 



-5/2 



(4.20) 
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As we see the post- Newtonian order of magnitude of AT is c +5 , that is the inverse of the order c~ 5 of radiation 
reaction effects. As a matter of fact, Af is a large number, approximately equal to 1.6 x 10 4 in the case of two neutron 
stars between 10 and 1000 Hz (roughly the frequency bandwidth of the detector VIRGO). Data analysts of detectors 
have estimated that, in order not to suffer a too severe reduction of signal to noise, one should monitor the phase 
evolution with an accuracy comparable to one gravitational-wave cycle (i.e. SAT ~ 1) o r bet ter. Now it is clear, from 
a post-Newtonian point of view, that since the 'Newtonian' number of cycles given by ( 4.2C ) is formally of order c +5 , 
any post-Newtonian correction therein which is larger than order c -5 is expected to contribute to the phase evolution 
more than that allowed by the previous estimate. Therefore, one expects that in order to construct accurate templates 
it will be necessary to include into the phase the post-Newtonian corrections up to at least the 2.5PN or 1/c 5 order. 
This expectation has been confirmed by various studies |2l|-|2~4|| which showed that in advanced detectors the 2.5PN 
or, better, the 3PN approximation is required in the case of inspiralling neutron star binaries. Notice that 3PN here 
means 3PN in the centre-of-mass energy E, which is deduced from the 3PN equations of motion, as well as in the total 
flux C, which is computed from a 3PN wave-generation formalism. For the moment the phase has been completed to 
the 2.5PN order ||-|^Jll; the 3PN order is still incomplete (but, see OM|9||). 



V. POST-NEWTONIAN WAVE GENERATION 



A. Field equations 

We consider a general compact-support stress-energy tensor T^ v describing the isolated source, and we look for 
the solutions, in the form of a (formal) post-Newtonian expansion, of the Einstein field equations, 

RM u _ I v R = ^—T^, (5.1) 
2 y c 4 K ' 

and thus also of their consequence, the equations of motion V „T Ml/ — of the source. We impose the condition of 
harmonic coordinates, i.e. the gauge condition 

d v h? v = 0; hf™ = V^gg^ - r? M ", (5.2) 

where g and g^ denote the determinant and inverse of the covariant metric g^, and where if LV is a Minkowski metric: 



jf- v = diag(— 1, 1, 1, 1). Then the Einstein field equations (5.1) can be replaced by the so-called relaxed equations, 
which take the form of simple wave equations, 

ahlxv = HtoG 
c 4 

where the box operator is the flat d'Alembertian □ = rj^d^dv, and where the source term t^ v can be viewed as the 
stress energy pseudotensor of the matter and gravitational fields in harmonic coordinates. It is given by 

r ^ = | g | T ^ + _|L A /- (5.4) 

is not a generally-covariant tensor, but only a Lorentz tensor relative to the Minkowski metric fy M „. As a 
consequence of the gauge condition ([T^), t^ v is conserved in the usual sense, 

d v T» v = (5.5) 

(this is equivalent to V ' JT^ V = 0). The gravitational source term K^ v is a quite complicated, highly nonlinear 
(quadratic at least) functional of h y,v and its first- and s eco nd-spa cetim c derivatives. 



We supplement the resolution of the field equations (5.2) and (|5.3[) by the requirement that the source does not 



receive any radiation from other sources located very far away. Such a requirement of 'no-incoming radiation' is to 
be imposed at Minkowskian past null infinity (taking advantage of the presence of the Minkowski metric rj^); this 
corresponds to the limit r = |x| — > +oo with t + r/c = constant. (Please do not confuse this r with the same r 
denoting the separation between the two bodies in section |iy|.) The precise formulation of the no- incoming radiation 
condition is 



lim 

t-F — =conatant 



(x,f)=0. (5.6) 
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In addition, rd\h^ v should be bounded in the same limit. Actually we often adopt, for technical reasons, the more 
restrictive condition that the field is stati onar y before some finite instant — T in the past (refer to || for details). 
With the no-incoming radiation condition ( |5.5| ) or (E^H) we transform the differential Einstein equation ( |5.3| ) into the 
equivalent integro-differential system 

hT = ^□hV", (5.7) 
where D^ 1 denotes the standard retarded inverse d'Alembertian given by 

(□-V)(x,t)=-i- J ^_ r (x',t-|x-x'|/c). (5.8) 



B. Source moments 



In this section we shall solve the field equations ( p.2| ) and (5.3) in the exterior of the isolated source by means 
of a multipole expansion, parametrized by some appropriate source multipole moments. The particularity of the 
moments we shall obtain, is that they are defined from the formal post-Newtonian expansion of the pseudotensor 
t^ v , supposing that the latter expansion can be iterated to any order. Therefore, these source multipole moments 
are physically valid only in the case of a slowly- moving source (slow internal velocities; weak stresses). The general 
structure of the post-Newtonian expansion involves besides the usual powers of 1/c some arbitrary powers of the 
logarithm of c, say 

r^(t, x, c) = ^^(*- x )> ( 5 - 9 ) 

where the overbar denotes the formal post-Newtonian expansion, and where are the functional coefficients of the 
expansion (p, q are integers, including zero). Now, the general multipole expansion of the metric field ft'* 1 ', denoted 
by M^W), is found by requiring that when re-developed into the near-zone, i.e. in the limit where r/c — > (this is 
equivalent with the formal re-expansion when c — > oo), it matches with the multipole expansion of the post-Newtonian 
expansion h (whose structure is similar to ([[1])) in the sense of the mathematical technics of matched asymptotic 
expansions. We find 00 that the multipole expansion M^h^) satisfying the matching is uniquely determined, 
and is composed of the sum of two terms, 

M(hn = finite part □ R 1 [M(A'-)] - -J E V"^ r^ (t ~ T,c) ' (5 ' 10) 

The first term, in which D^ 1 is the flat retarded operator (5.8), is a particular solution of the Einstein field equations 
in vacuum (outside the source), i.e. it satisfies Uh^ t = A4(A^"). The second term is a retarded solution of the 
source-free homogeneous wave equation, i.e. ^ft^cm = 0- We denote 8l — . . . di t where L = u ... is a multi- 



index composed of / indices; the I summations over the indices are not indicated in ( 5.1C ). The 'multipole 

moments'parametrizing this homogeneous solution are given explicitly by (with u = t — r/c) 

Tl v (u) = finite party d 3 xx L J dz ^(z)r^(x, u + z\x\/c), (5.11) 

where the inte gran d con tains the post- Newtonian expansion of the pseudostress-energy tensor r'" y , whose structure 



reads like (5.S). In ( 5.11 ), we denote the symmetric-trace- free (STF) projection of the product of I vectors x l with a 



hat, so that xl = STF(x L ), with x L = x 11 . . . x 11 and L = i\ . . . if, for instance, Xij — XiXj — 



8i(z) is given by 



and satisfies the properties 



S l (z)= { ^±^(l-zy, (5.12) 
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dz5i(z) — 1; lira 6i{z) — S(z) 



l — >-\-oo 



(5.13) 



(where 5(z) is the Dirac measure). Both terms in (5.10) involve an operation of taking a finite part. This finite part 
can be defined precisely by means of an analytic continuation (sec Jbij for details), but it is in fact basically equivalent 
to taking the finite part of a divergent integral in the sense of Hadamard [jl8| . Notice, in particular, that the finite 
part in the expression of the multipole moments ( [5. 11 ) deals with the behaviour of the integral at infinity: r — > oo 
(without the finite part the integral would be divergent because of the factor xl — r l riL in the integrand and the fact 
that the pseu doten so r T^" is not of compact support). 

The result (5.1C)-( |5.11 ) permits us to define a very convenient notion of the source multipole moments (by oppo- 
sition to the radiative moments defined below). Quite naturally, the source moments are constructed from the ten 
components of the tensorial function T^ v (u). Among these components four can be eliminated using the harmonic 
gauge condition (5.2), so in the end we find only six independent source multipole moments. Furthermore, it can be 
shown that by changing the harmonic gauge in the exterior zone one can further reduce the number of independent 
moments to only two. Here we shall report the result for the 'main' multipole moments of the source, which are the 
mass-type moment II and current-type Jl (the other moments play a small role starting only at highorder in the 
post-Newtonian expansion). We have fli3] 

I L {u) = finite part J d x J dz y>fi L E - c 2^ + ^(21 + $j S l+ 1 ^ iL ® t ^ i 



2(21 + 1) 



c 4 (Z + l)(Z + 2)(2Z + 5) 
Jl(u) = finite part J d 3 x J dz e ah<il j<Jj2C£_i> E 



21 + 1 



c 2 (Z + 2)(2Z + 3) 



<5/+iXL_i> ac 9 t S 



be 



(5.14) 



(5.15) 



Here the integrand is evaluated at the instant u + z\x\/c, e a b c is the Levi-Civita symbol, (L) is the STF projection, 
and we employ the notation 



2 ' i ^l] 

c? c 



(5.16) 



(with ' 



5ijT tJ ). The multipole moments II, Jl are valid formally up to any post-Newtonian order, and constitute 



a generalization in the nonlinear theory of the usual mass and current Newtonian moments (see, ]14j| for details). It 
can be checked that, when considered at the 1PN order, these moments agree with the different expressions obtained 
in j^] (case of mass moments) and in pc| ] (current moments). 



C. Radiative moments 



In linearized theory, where we can neglect the gravitational source term A MI/ in ( p.4| ) , as well as the first term in (5.10), 
the source multipole moments coincide with the so-called radiative multipole moments, defined as the coefficients of 
the multipole expansion of the 1/r term i n the distance to the source at retarded times t — r/c = constant. However 



in full nonlinear theory, the first term in (5.1C) will bring another contribution to the 1/r term at future null infinity. 
Therefore, the source multipole moments are not the 'measured' ones at infinity, and so they must be related to the 
real observables of the field at infinity which are constituted by the radiative moments. It has been known for a long 
time that the harmonic coordinates do not belong to the class of Bondi coordinate systems at infinity, because the 
expansion of the harmonic metric when r — > oo with t — r/c = constant involves, in addition to the normal powers 
of 1/r, some powers of the logarithm of r. Let us change the coordinates from harmonic to some Bondi-type or 
'radiative' coordinates (X, T) such that the metric admits a power-like expansion without logarithms when R — > oo 
with T — R/c = constant and R = |X| (it can be shown that the condition to be satisfied by the radiative coordinate 
system is that the retarded time T — R/c becomes asymptotically null at infinity). For the purpose of deriving the 



formula ( 5.20| ) below it is sufficient to transform the coordinates according to 



T — — = t -- - —j— to (— ) , (5.17) 

c c c a \ro 
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where M denotes the ADM mass of the source and ro is a gauge constant. In radiative coordinates it is easy to 
decompose the 1/ R term of the metric into multipoles and to define in that way the radiative multipole moments Ul 
(mass-type; where L — i\ . . .it with I > 2) and Vl (current- type; with I > 2). (Actually, it is often simpler to bypass 
the need for transforming the coordinates from harmonic to radiative by considering directly the TT projection of the 
spatial components of the harmonic metric at infinity.) The formula for the definition of the radiative moments is 



,TT 



AG 



i=i 



N L -2U abL - 2 {T - R/c) 



c(l 



^N cL ^e cd{a V b)dL _ 2 (T - R/c) I + O 



(5.18) 



where N is the vector Ni = N' 1 = X l /R (for instance Nl~2 = . . . Ni t _ 2 ), and Vij ab denotes the TT projector 



Vijab = (S ia - NiN a )(6 jb - NjNb) - - NiNjXSeb - N a N b ). 



(5.19) 



In the limit of linearized gravity the radiative multipole moments Ul, Vl agree with the Zth time derivatives of the 
source moments II, Jl- Let us give, without proof, the result for the expression of the radiative mass-quadrupolc 
moment [Zy including relativistic corrections up to the 3PN or 1/c 6 order inclusively p^Jl^ |. The calculation involves 
implementing explicitly a post-Minkowskian algorithm defined in || for the computation of the nonlinearities due to 
the first term of (Ol. We find (U = T - R/c) 



U ij {U) = M ( >?{U) + 2 



GM 



G 



dvM^(U-v) 



dv [M<%M$ a ]{U-v 

(5) ^- a(U) 



2r 



o 

(i) 



-Mi%MY>(U) + -M" 



-M {3) M {2) 
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12 

] a(U) 



<i M 3>° 



-s ab<i M^ a J b (U) 



dv M (5) 



(u- 



cv 
2n> 



57 
70 



In 



cv 
2^ 



124 627 
44100 



(5.20) 



The superscript (n) denotes n time derivations. The quadrupole moment My entering this formula is closely related 
to the source quadrupole , 



Mi. 



(5.21) 



where K is the Newtonian moment of inertia (see equation (4.24) injEjl; we are using here a mass-centred frame so 
that the mass-dipole moment Ii is zero). The Newtonian term in ( |5.20| ) corresponds to the quadrupole formalism. 
Next, there is a quadratic nonlinear correction term with multipole interaction M x M^- which represents the effect 
of tails of gravitational waves (scattering of linear waves off the spacetime curvature generated by the mass M). This 
correction is of order 1/c 3 or 1.5PN and takes the form of a non-local integral with logarithmic kernel p0[ . It is 
responsible notably for the term proportional to 7ft 1 / 4 in the formula for the phase ( 3.13 ) below. The next correction, 
of order 1/c 5 or 2.5PN, is constituted by quadratic interactions between two mass-quadrupoles, and between a mass- 
quadrupole and the constant current dipole Jl2| . This term contains also a non-local integral, which is due to the 
radiation of g ravita tional waves by the distribution of the stress-energy of linear waves plp^ , |30| , p^| . Finally, at the 
3PN order in ( [5.20 ) the first cubic nonlinear interaction appears, which is of the type (M x M x M^Jj and corresponds 
to the tails generated by the tails themselves [^31 . 



VI. INSPIRAL BINARY WAVEFORM 



To conclude, let us give (without proof) the result for the two polarization waveforms h+(t) and h x (t) of the 
inspiralling compact binary developed to 2PN order in the amplitude and to 2.5PN order in the phase. The calculation 
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was done by Blanchet, Damour, Iyer, Will and Wiseman |25| , |26| , [T5| , p7| , ^3| , based on the formalism reviewed in section |v| 
and, independently, on that defined in [filif . Following |33| we present the polarization waveforms in a form which is 
ready for use in the data analysis of binary inspirals in the detectors VIRGO and LIGO (the analysis will be based 
on the optimal filtering technique reviewed in section |ll|). We find, extending the Newtonian formulae in section 

_ 2G> ( Gmu \^ ,] 



III. 



x{< x +- 1/2 <x 2) +x 3 M% 2) +^l 2) x }, (6.1) 
where the various post-Newtonian terms, ordered by x, are given for the plus polarization by 

iz| 0) = -(l + c 2 )cos2V>, (6.2) 

H +' 2) =~j ^[(5 + c 2 )cosV>-9(l + c 2 )cos3^], (6.3) 

H { + ] = ±[(19 + 9c? - 2c?) - i/(19 - 11c? - 6c?)] cos 2^ 

-§.s 2 (l + c 2 )(l-3^)cos4V>, (6.4) 

-^[(73 + 40c 2 - 9c?) - 2^(25 - 8c 2 - 9c?)] cos 3-0 

+2|5(1 - 2^)s 2 (l + c 2 ) cos5?/>} - 2tt(1 + c\) cos2V>, (6.5) 

H + ] = m t( 22 + 396c2 + 145c ^ - 5c ^) + f^(706 - 216c 2 - 251^ + 15c?) 
-5z/ 2 (98 - 108c 2 + 7cf + 5ef)] cos 2^ 
+ prs?[(59 + 35c 2 - 8c?) - |^(131 + 59c 2 - 24c?) 



3' 

+5i/ 2 (21 - 3c 2 - 8c?)] cos 4^ 
-§(1 - hv + 5i/ 2 )s? (1 + c 2 ) cos6?A 

+ ^- — {[11 + 7c 2 + 10(5 + c 2 ) ln2] sin7/> - 5tt(5 + c 2 ) cos?A 
40 m 

-27[7- 10 ln(3/2)](l + c 2 )sin3V> + 135tt(1 + c 2 )cos3Vj}, 



and for the cross-polarization by 



# 5 = | [(17 - 4c 2 ) - K13 - 12c 2 )] sin2^ - |(1 - 3^)c, S 2 t 



^(3/2) = £i£i^ {[(63 _ 5c 2 } _ 2jy(2 3 _ 5c2)] sin ^ 
96 to 



ff(, 2) = |i[(68 + 226c- - 15c?) + §i/(572 - 490c 2 + 45c 4 ^ 

7nJ i i r„4\ 



60 1 

-5^ z (56 - 70c- + 15c?)] sin 2^ 
+^c 4 s 2 [(55 - 12c 2 ) - §1^(119 - 36c 2 ) + 5^ 2 (17 - 12c 2 )] sin4V> 



(6.6) 



H { ° } = -2 Cl sin2^, (6.7) 

H^y 2 " 1 = — -SjCj^^ [sin i/) — 3 sin 3^], (6.8) 
4 m 



(6.9) 



-f [(67- 15c 2 ) - 2^(19- 15c 2 )] sin3?A 

+^(1 - 2v)s\ sin 5^} - 4tt Ci sin 2ip, (6.10) 



§(l-5z/ + 5i/ 2 )ciS? sin 6^ 

3 $7Tl 

—SiCi {[3 -I- 10 In 2] cos?/; + 5-Ksmijj 

20 m 

9[7- 101n(3/2)] cos3?/> - 457rsin3V>}. (6.11) 
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The n otation is consistent with sections III and IV. In particular, the post-Newtonian parameter x is defined by 
( 4.13 ). We use the shorthands Cj = cosi and Sj = sini where i is the inclination angle. The basic phase variable ip 
entering the waveforms is defined by 



2Gmuj ( lu 
ip = 4> 5 — In — 

C J \CJQ 



(6.12) 



where 4> is the actual orbital phase of the binary, and where ujq can be chosen as the seismic cut-off of the detector (see 
|33| fo r details). As for the phase evolution </>(t), it is given up to 2.5PN order, generalizing the Newtonian formula 
MM, by 



4>{t) = 0o 



T 5/8 



1 

V 

9 275 495 
14 450 688 
38 645 



3715 55 



8064 96 
284 875 



r 3/8 



-1/4 



1855 



258 048 
15 



2048 



-1/8 



172 032 2048 



V 7T In T 



(6.13) 



where <f>o is a constant and where we recall that the dimensionless time variable r was given by (4.17). The frequency 
is equal to the time derivative of (6.13), hence 



w(t) 



8Gm 



.-3/8 



743 11 



2688 32 
56 975 



1855 099 



14 450 688 258 048 



371 

2048 1 
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.-7/8 



7729 3 
21504 ~ 256 Z 



(6.14) 



We have checked that both waveforms ([T^-Q and phase/frequency (|6Tll|)-(p4|) agree in the test mass limit 
v — ► with the results of linear black hole perturbations as given by Tagoshi and Sasaki pS . 
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